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ABSTRACT 

Wavelet  theory  provides  an  attractive  approach  to  signal  and  image  compression.  This  work  investigates  a 
new  approach  for  wavelet  transform  coefficient  selection  for  efficient  image  compression.  For  a  desired 
image  compression  ratio  (50:1),  wavelet  scale  thresholds  are  derived  via  a  multiagent  stochastic 
optimization  process.  Previous  work  has  demonstrated  an  interscale  relationship  between  the  stochastically 
optimized  wavelet  coefficient  thresholds.  Based  on  the  experimental  results,  a  deterministic  wavelet 
coefficient  selection  criterion  is  hypothesized  and  the  constants  of  the  equation  statistically  derived. 
Experimental  results  of  the  stochastic  optimization  and  deterministic  approaches  are  compared  and 
contrasted  with  results  from  previously  published  wavelet  coefficient  threshold  strategies. 

Keywords:  wavelet  transform,  image  compression,  stochastic  optimization. 

1.  INTRODUCTION 

Wavelets  originated  from  the  analysis  of  seismic  data  by  J.  Morlet.^  A  wavelet  is  a  function  y/  whose 
dilations  and  translations 


=  a>Q,  (1) 

form  a  basis  in  The  continuous  wavelet  transform  maps  a  function  /  to  a  set  of  wavelet  coefficients 

which  vary  continuously  over  {a,b) 

^  J  m  (2) 

Given  an  orthonormal  wavelet  basis,  any  Is  function  /  can  be  represented  as  a  (possibly  infinite) 
summation 

=  (3) 


where  the  wavelet  coefficients  are  defined  as  the  scalar  inner  product 


<,  =  ]m¥St)dt- 


(4) 


The  wavelet  coefficients  correspond  to  the  projection  of  f{t)  onto  the  coefficient’s  associated  wavelet. 
Wavelet  basis  sets  can  be  orthogonal  or  non-orthogonal,  and  each  wavelet  can  have  compact  or  infinite 
support.  For  image  compression,  one  does  not  want  to  transform  the  image  into  more  coefficients  than 
necessary  to  completely  characterize  the  image.  Orthogonal  wavelets,  such  as  the  Haar  and  Daubechies 
wavelets,  provide  an  efficient  and  economical  approach  to  the  decomposing  of  images  into 
multiresolutional  representations."^  Therefore,  an  orthogonal  wavelet  (the  four  coefficient  Daubechies 
wavelet)  is  used  in  this  research. 


Traditionally,  image  transform  coding  or  compression  methods,  such  as  the  discrete  cosine  transform,  rely 
on  high  energy  compaction  with  an  m  x  n  passband  coefficient  selection  criterion.  The  wavelet  tra^isform 
generally  results  in  poor  energy  compaction,  so  an  alternative  strategy  is  required  for  wavelet  coefficient 
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selection.  This  paper  investigates  a  multiagent  stochastic  optimization  approach  for  wavelet  coefficient 
selection,  and  a  statistical  derivation  of  a  deterministic  equation  for  wavelet  coefficient  selection.  This 
derivation  is  based  on  empirical  data  derived  from  the  stochastic  optimization  of  the  coefficient  thresholds. 

LI.  Multiresolution  coding  theory 

Performing  the  Fourier  transform  decomposes  an  image  into  a  frequency  representation  and  allows  an 
analyst  to  determine  if  a  particular  frequency  is  present,  but  does  not  provide  any  spatial  information.  The 
wavelet  transform,  unlike  the  Fourier  transform,  produces  a  spatial-frequency  representation  that  has  good 
localization  in  both  domains. 

The  wavelet  transform  filters  the  data  and  results  in  a  lowpass,  or  blurred,  version  of  the  original,  and  a 
highpass,  or  detailed,  version.  The  lowpass  image  is  recursively  filtered  to  produce  a  blurred  and  detail 
image  of  itself  Because  the  wavelet  transform  used  in  this  work  is  separable,  the  image  wavelet  transform 
can  be  implemented  as  successive  one-dimensional  transforms  of  the  row  and  columns.  The  resulting 
image  has  a  horizontal  detailed  image,  a  vertical  detailed  image,  a  diagonal  detailed  image,  and  a  blurred 
image.  Each  transform  consists  of  a  wavelet  dilation,  or  wavelet  resolution. 

We  are  interested  in  a  discrete  sublattice  of  the  continuous  wavelet  since  image  data  is  discrete  in  nature. 
Thus,  the  discrete  wavelet  is 


^  x-nb^a;^ 


(5) 


where  and  are  constants.  In  our  case,  =  2  and  ^  •  Rewriting  the  wavelet  transform  in  terms 
of  the  inner  products  of  the  finer  resolution  wavelet  transform  shows  that  the  discrete  multiresolutional 
decomposition  can  be  performed  using  the  filters  h  and  g,  given  that 

H{0)  =  1 

\H{cof+\H{co^Kf=\  (6) 

where  H{co)  represents  the  Fourier  transform  of  the  filter  The  corresponding  highpass  filter  g  is 

defined  via 


G{(o)  =  e-^^H\o:)  +  7i)  (7) 

where  //*(•)  is  the  complex  conjugate  of  //(•)>  and  G(a))  is  the  Fourier  transform  of  g. 

Directly  implementing  the  wavelet  transform  as  given  in  Eq.  2  is  inefficient.  A  multiresolutional 
decomposition  algorithm  has  been  developed  that  accomplishes  the  wavelet  transform  in  order 
log  A  operations  where  N  is  the  length  of  the  data.*^  The  algorithm  introduces  an  auxiliary  function 
(p{x)  related  to  Y{x)  by 

^fix)  =  ^h{k)(l>{x-k)  (8) 

k 

where  (l>(x)  satisfies 

(p{x)  =  Y,8{k)<P(2x-k).  (9) 


1.2.  Wavelet  compression  techniques 

Wavelet  image  coding  techniques  include  wavelet^  and  coefficient^  selection  approaches,  as  well  as  hybrid 
schemes  which  combine  wavelets  with  other  signal  representations. In  a  coefficient  selection  approach, 
image  compression  is  achieved  by  retaining  a  subset  of  coefficients  from  the  wavelet  representation,  and 
zeroing  the  remaining  coefficients.  The  number  of  coefficients  selected  is  based  on  the  maximum  mean- 
squared  error  (MSE)  acceptable,  the  compression  ratio  required,  or  a  combination  of  both.  The  coefficients 
selected  are  then  coded  using  an  entropy-based  encoder  that  stores  the  value  of  the  coefficient  ,and  its 
location  within  the  current  scale. 


Different  sets  of  retained  coefficients  can  result  in  drastically  different  image  reconstructions.  The 
traditional  method  for  selection  is  simply  a  passband  filter  of  some  height  and  width  that  passes  a 
rectangular  region  of  the  transformed  data  and  sets  the  rest  to  zero.  DeVore  et  al.  show  that  the  wavelet 
transform  does  not  compact  energy  into  a  particular  region,  making  the  passband  approach  suboptimal. 
DeVore  presents  a  wavelet  coefficient  selection  algorithm  that  keeps  the  highest  N  coefficients  and 
discards  the  rest.  Argast  et  al.  demonstrate  that  coefficients  from  lower  resolution  scales  affect  a  larger 
part  of  the  image  than  coefficients  in  the  high  resolution  scales,  Argast  presents  an  algorithm  for  selecting 
coefficients  based  on  a  logarithmic  threshold,  where  the  threshold  of  each  scale  is  twice  the  threshold  of  the 
previous  lower  resolution  scale.  This  compression  strategy  of  Argast  is  given  as: 


for  J  =  S  to  \  [ 

for  each  element  k  in  scale  J 
if  CR%-  GOAL  then  stop 
else 


Jo 

set  c,  -  < 

}  end 


if  K|<A-2^ 
otherwise 


where  c^  is  a  transform  coefficient,  CR%  is  the  current  compression  percentage,  is  a  constant  weight 
factor,  and  J  and  S  are  respectively  the  current  and  total  wavelet  scales.  The  constant  A  is  set  (initially 
set  to  0.01  in  their  paper)  and  then  iteratively  adjusted  so  that  the  compression  ratio  can  be  met.  For  a 
selected  compression  ratio,  this  approach  results  in  a  lower  mean-squared  error  in  the  reconstructed  images 
than  DeVore’s  algorithm. 

Previous  preliminary  work^  has  demonstrated  an  interscale  relationship  between  an  image’s  stochastically 
derived  wavelet  coefficient  thresholds.  This  paper  increases  the  number  of  images  stochastically 
optimized,  and  for  a  selected  compression  ratio,  statistically  derives  from  empirical  results  a  deterministic 
approach  for  wavelet  coefficient  selection. 

2.  STOCHASTIC  OPTIMIZATION 

Since  the  1950’s,  random  (or  stochastic)  search  techniques  have  been  used  for  function  optimization. 
Random  search  strategies  are  competitive  with  traditional  search  strategies  (such  as  gradient  search 
techniques)  when  the  cost  or  objective  function  is  expensive  or  difficult  to  compute,  or  when  the  function  to 
be  minimized  has  many  suboptimal  solutions  (local  minima).  Other  advantages  include  the  ease  of 
programming,  inexpensive  realization  of  possible  solutions,  as  well  as  flexibility  in  the  expression  of  the 
criterion  function. 

Stochastic  optimization  techniques  are  based  on  either  a  single  point  or  multiple  agent  algorithms.  Single 
point  algorithms  include  the  random  walk,  the  creeping  random  method, and  the  method  of  Solis  and 
Wets.^^  Multiple  agent  stochastic  search  algorithms,  such  as  genetic  algorithms,  evolutionary  strategies, 
and  evolutionary  programming  (EP),  are  becoming  well  known  for  their  optimization  properties.  A 
thorough  discussion  of  the  EP  algorithm,  and  other  evolutionary  algorithms,  is  given  by  Back  and 
Schwefel. 

In  EP,  a  population  of  models  (solutions)  generate  new  models  via  a  mutation  process.  The  new  models 
compete  with  the  original  (parent)  population  for  survival  to  the  next  iteration.  The  mutation  process  is  a 
random  perturbation  of  the  parameters,  with  the  magnitude  of  the  perturbation  generally  tied  to  the  fitness 
of  the  parent.  Although  useful  in  escaping  local  minima,  the  random  nature  of  the  stochastic  mutation 
process  can  be  inefficient.  To  increase  convergence  efficiency,  variants  of  the  EP  algorithm  have  been 
developed.  ^^’20  The  presented  work  imbeds  the  method  of  Solis  and  Wets  into  the  traditional  evolutionary 
programming  paradigm.  The  parent  models  not  only  produce  a  single  offspring  by  way  of  the  mutation 
process,  but  are  also  modified  via  the  Solis  and  Wets  algorithm.  Figure  1  illustrates  this  hybrid 
optimization  approach. 


component  corresponds  to  the  final  low-pass  (DC)  component  of  the  wavelet  transform).  Figure  2  displays 
a  multiresolutional  wavelet  transform,  and  the  threshold  array  indices  corresponding  to  the  wavelet 
transform  scales. 


Figure  2.  Multiresolutional  wavelet  transform  and  associated  scales. 

For  each  image,  the  compression  ratio  was  set  to  50:1  (CA%  =  2.0%).  The  four  coefficient  Daubechies 
wavelet  is  used  for  the  compression  process.  Since  we  are  currently  interested  in  the  coefficient  thresholds, 
optimization  of  the  wavelet  is  beyond  the  scope  of  this  paper.  Figure  3  displays  the  MSE  and  percentage  of 
selected  coefficients  of  the  best  coefficient  threshold  strategy  for  two  images  during  the  stochastic 
optimization  process.  As  shown  in  the  figures,  the  best  solution  of  the  evolutionary  programming  process 
quickly  achieves  the  desired  compression  ratio. 


Iteration  Iteration 


MSE 

Compression 


Figure  3.  Stochastic  optimization  performance  for  two  images. 

For  all  of  the  images  compressed,  the  coefficient  threshold  values  of  the  best  compression  strategy  (the  best 
T  vector)  derived  after  600  iterations  are  tabulated  in  Table  1,  The  stochastically  derived  thresholds  for 
several  of  these  images  are  also  shown  in  Figure  4.  It  is  important  to  note  that  for  all  images,  the  values  for 
the  lower  scales  are  below  all  wavelet  coefficient  values  in  those  scales,  and  therefore  all  coefficients  in 
these  scales  are  kept  for  image  reconstruction  by  the  compression  strategy.  Table  1  italicizes  the  threshold 
values  which  are  less  than  the  minimum  wavelet  coefficient  value  in  the  respective  scale.  As  shown  in  the 
table  and  figure,  the  higher  scale  thresholds  generally  demonstrate  an  interscale  relationship,  with  the  value 
of  the  higher  scale  roughly  double  the  value  of  the  lower. 


Table  1.  Derived  wavelet  coefficient  thresholds  (italicized  thresholds  <  min,  scale  coefficient  value). 
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Figure  4.  Wavelet  coefficient  thresholds  for  several  images. 


To  measure  if  the  coefficient  threshold  values  derived  from  the  optimization  process  are  related  to  the 
variance  or  power*  of  an  image,  the  intensity  histogram  of  the  "watch”  image  was  manipulated  and  the 
resulting  new  images  where  stored.  The  coefficient  thresholds  for  each  version  were  stochastically 
optimized,  and  the  resulting  thresholds  are  displayed  in  Table  2  and  Figure  5.  The  amount  of  energy  and 
power  of  each  image  were  also  calculated,  and  each  image’s  associated  power  is  used  to  label  the  images. 
The  figure  and  table  demonstrate  that  the  thresholds  are  functionally  related  (in  some  way)  to  the  energy  or 
power  of  the  image,  with  images  with  a  higher  power  content  requiring  higher  coefficient  thresholds  for 
efficient  image  compression. 


Table  2.  Thresholds  for  watch  image  versions  {italicized  thresholds  <  min,  coefficient  value). 
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*  The  variance  or  power  of  an  image  is  defined  as  ,  where  /(/y)is  the  pixel  intensity  aK 

^ j  ij 

location  (/,;)  and  p.,  is  the  mean  intensity  value  of  the  image. 


Figure  5.  Coefficient  thresholds  for  an  image  at  various  levels  of  power. 

5.  STATISTICAL  ANALYSIS  OF  HYPOTHESIZED  COMPRESSION  SCHEME 

In  this  section,  we  hypothesize  the  nature  of  the  relationship  suggested  by  the  results  of  the  multiagent 
stochastic  optimization  process,  and  derive  a  deterministic  equation  based  on  the  statistics  derived  from  the 
data.  For  a  given  compression  ratio,  we  propose  a  coefficient  selection  criterion  of  the  form 

T=K(CR)-fia^)-2'  (12) 

where  T.  is  again  the  coefficient  threshold  associated  with  the  ith  wavelet  scale,  K{CR)  is  a  function  of 
the  compression  ratio,  and  /(O'*)  is  a  function  of  the  power  content  of  the  image.  In  this  work,  we 
hypothesize  the  functional  form  of  f{G^)  to  be 


(13) 

where  B  is  a  constant.  In  this  investigation,  a  single  compression  ratio  has  been  applied  to  all  images,  and 
therefore  K  is  considered  a  constant.  The  hypothesized  equation  can  be  written  as 

T.,=K{a~yT  (14) 

We  will  now  apply  regression  analysis  techniques  to  calculate  a  value  for  B  based  on  the  thresholds 
derived  in  the  stochastic  optimization  process.  For  both  constants  B  and  K ,  only  the  coefficient  thresholds 
from  scales  four  and  above,  which  actually  zero  coefficients  in  the  scale,  are  used  for  the  statistical 
derivation.  Threshold  values  in  scales  zero  through  three  are  not  used  in  the  analysis,  since  any  threshold 
value  below  the  minimum  coefficient  value  of  a  scale  (i.e.  any  value  T.  such  that  0  <  7.  <min{lc^J}) 

results  in  the  same  image  reconstruction. 

5.1.  Determination  of  the  relationship  between  threshold  value  and  power 

The  intensity-adjusted  "watch"  images,  and  the  stochastically  derived  coefficient  thresholds,  are  used  to 
determine  the  relationship  between  image  power  and  the  threshold  values.  Dividing  both  sides  of  Eq.  14  by 
2'  and  taking  the  logarithm  of  both  sides  gives 

log(f)  =  51og(6'^)  +  logA:.  (15) 

This  is  a  linear  equation  y  =  5jc-hc,  with  y  =  log(y),  .r  =  log(G'“),  and  c  =  log  A',  and  linear  regression 

techniques'^  are  applied  for  the  estimate  of  B.  Given  the  assumption  that  the  are  samples  are  independent 
and  normally  distributed  about  the  true  values  also  allows  for  statistical  inference  and  validation  tests  to  be 
performed.  The  maximum  likelihood  regression  estimate  of  B  is  given  as 


Figure  6  displays  the  sample  pairs  y  =  log(y)  vs.  x  =  log(P),  and  line  determined  by  the  linear  regression 
estimates  for  the  threshold  coefficients  under  analysis. 


Figure  6.  Relationship  between  and  y  =  log(|-)  and  x  =  logCa^). 

The  maximum  likelihood  value  for  the  slope  is  computed  to  be  B  =  0.4714,  with  an  associated  coefficient 

of  determination  for  the  regression  line  is  equal  to  =91.37%.  The  value  of  B  is  very  close  to  0.5,  and 
assuming  the  normal  model  allows  the  use  of  the  F-statistic  associated  with  the  hypothesis  that 
vs.  H^'.B^Yi.  Computing  an  a  =  0. 10  level  of  significance  test,  the  corresponding  F-statistic  requires  one 
to  reject  H^:B-Yi  if  the  statistic  is  greater  than  Fy  (1, 40)  =  2.84.  The  F-statistic  is  given  as 


(0.5-0.4714)' 


0. 


=  1.560 


(17) 
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which  is  less  than  (1, 40) .  So  one  fails  to  reject  that  H^\B  =  Yi  based  on  the  data  analyzed.  This  reduces 
the  functional  from  of  Eq.  13  to 


(18) 


Since  5  =  0.5  provides  a  simple  interpretation,  with  coefficient  threshold  values  directly  related  to  the 
standard  deviation  of  the  image  intensities,  we  will  use  5  =  0.5  for  the  determination  of  the  constant  K . 

5.2.  Determination  of  the  constant  K 

Given  the  derived  value  of  B,  our  hypothesized  equation  (Eq.  14)  is  now 

r.  =  KaT  (19) 

Dividing  both  sides  by  a2'  derives  a  value  of  K.  for  each  image  and  scale  under  analysis,  and  these  values 
are  used  as  a  sample  to  derive  an  estimate  of  K .  Nine  images  of  Table  1  are  used  for  the  computation  of 
K ,  with  the  "Lenna”  image  reserved  for  algorithm  testing  and  comparison.  It  was  found  that  the  mean 
value  of  K  did  not  achieve  adequate  compression  on  the  test  image,  since  the  mean  will  “on  average" 
produce  the  desired  compression.  To  more  adequately  achieve  the  desired  compression  ratio,  one-half  of  a 
standard  deviation  of  the  {/C.}  sample  is  added  to  the  mean  K  to  produce  K .  K  is.  simply  computed  as 


K  =  K^^ 
2 


(20) 


Note  that  the  highest  scale  threshold  value  of  image  ”woman2”  is  a  wild  point,  and  this  value  was  not  used 
in  the  computation  of  K .  For  the  images  analyzed,  the  estimate  K  is  computed  dLS  K  =  9.360  •  10“*,  and 
the  standard  deviation  of  {ATJ  is  s,^  =  6.217  •10'^.  The  median  of  the  (A".)  sample,  a  more  robust  measure 

of  centrality  when  wild  points  exist,  is  in  close  agreement  with  A',with  =  9.604  •  10*^.  For  the  single 
compression  ratio  (50:1)  examined  in  this  paper,  our  hypothesized  deterministic  equation  is  finalized  as 

r,  =(1.247 •  10-' •o-)2'.  (21) 

Unlike  the  algorithms  of  DeVore  and  Argast  and  the  stochastically  optimization  approach,  the  coefficients 
of  Eq.  21  have  been  statistically  determined,  and  the  resulting  compression  ratio  for  an  image  using  Eq.  21 
will  not  be  exactly  50:1.  However,  this  equation  requires  neither  coefficient  sorting  nor  iterative 
optimization,  and  is  therefore  computationally  more  efficient  than  the  other  approaches.  The  next  section 
compares  this  deterministic  equation  with  the  algorithms  previously  discussed  in  Section  2  on  the  ’’Lenna" 
image. 


6.  COMPARISON  OF  IMAGE  COMPRESSION  APPROACHES 

The  "Lenna"  image,  not  used  for  the  statistical  determination  of  Eq.  21,  is  used  for  a  comparative  analysis 
of  the  wavelet  compression  algorithms.  Table  3  compares  the  algorithms  of  DeVore,  Argast,  the 
stochastically  determined  thresholds,  and  the  wavelet  coefficient  thresholds  calculated  using  the 
deterministic  equation  (Eq.  21),  on  the  "Lenna"  image. 


Table  3.  Comparison  of  algorithms  on  the  ’’Lenna"  image. 


esssbshhbi 

MSE 

%  of  Distortion 

DeVore 

135.24 

5.91  % 

12.29  dB 

50:1 

Argast 

75.80 

3.31  % 

14.80  dB 

50:1 

Stochastic 

68.26 

3.02  % 

15.19dB 

50:1 

Eq.  21 

75.87 

3.31  % 

14.80  dB 

56:1 

Unfortunately,  the  standard  deviation  of  the  sample  [K.}  is  on  the  same  order  of  magnitude  as  the  mean  K , 
and  this  indicates  that  the  estimates  of  K  vary  dramatically  in  the  sample.  This  variability  is  illustrated  in 
Table  4,  where  the  MSE  and  achieved  compression  ratio  are  compared  for  four  images.  As  shown  in  the 
table,  the  compression  ratio  for  the  equation  ranges  from  33:1  to  86:1.  This  is  due  to  the  large  deviation  of 
the  estimates  for  K  in  the  images  analyzed.  This  variability  severely  limits  the  usefulness  of  the 
deterministic  equation  when  achievement  of  a  compression  ratio  is  required,  and  brings  the  assumption  that 
A'  is  a  function  of  compression  ratio  into  question. 


Table  4.  Comparison  of  algorithm  MSE  and  compression  ratio  on  four  images. 


Image 

light 

bird 

fruits 

truck 

MSE 

CR 

MSE 

CR 

MSE 

CR 

MSE 

CR 

DeVore 

193.91 

50:1 

30.73 

50:1 

141.87 

50:1 

124.04 

50:1 

Argast 

129.38 

50:1 

19.81 

50:1 

80.13 

50:1 

80.74 

50:1 

Stochastic 

118.15 

50:1 

19.50 

50:1 

68.97 

50:1 

75.40 

50:1 

Eq.21 

152.28 

66:1 

34.30 

86:1 

83.68 

64:1 

58.26 

33:1 

As  shown  by  Tables  3  and  4,  the  wavelet  coefficient  thresholds  found  by  the  stochastic  optimization 
process  are  superior  to  all  other  approaches  discussed  in  this  work.  For  every  image  tested,  the  multiagent 
stochastic  process  produced  the  coefficient  threshold  strategy  with  the  minimum  MSE.  The  method  of 
Argast,  although  outperformed  by  our  stochastic  optimization  process,  was  significantly  better  than  the 
method  of  DeVore  for  all  images  tested. 


Figure  7  displays  the  original  and  restored  "Lenna"  images.  These  images  represent  the  results  from  the 
application  of  the  various  50:1  wavelet  compression  approaches,  each  using  the  four-coefficient  wavelet  of 
Daubechies. 


7.  CONCLUSIONS 

In  conclusion,  multiagent  stochastic  optimization  provides  a  valuable  tool  for  the  investigation  of  complex 
problems.  In  this  work,  we  have  developed  a  stochastic  optimization  technique  to  determine  the  wavelet 
coefficient  thresholds  which  best  compress  and  restore  an  image  for  a  given  compression  ratio.  With  the 
stochastically-derived  wavelet  coefficient  thresholds,  we  have  shown  evidence  of  a  strong  statistical 
relationship  between  threshold  magnitude  and  the  power  of  the  image.  For  every  image  analyzed  and 
tested,  the  compression  strategy  derived  via  our  stochastic  optimization  process  outperforms  the  iterative 
and  sorting  wavelet  coefficient  threshold  procedures,  producing  compression  strategies  that  achieve  the 
required  compression  with  the  consistently  minimal  MSE. 

A  deterministic  threshold  selection  equation,  based  on  the  stochastically-derived  thresholds,  was  derived. 
This  equation  has  shown  superior  capability  for  image  reconstructive  ability  when  compared  to  the 
approach  of  DeVore,  and  is  computationally  more  efficient  than  the  techniques  of  Argast  or  Waagen  et  al.^ 
However,  the  variability  of  the  hypothesized  constant  K  limits  the  equation’s  usefulness  for  image 
compression  to  situations  with  ”on-average”  compression  requirements. 

This  work  represents  ongoing  research,  and  future  work  will  investigate  the  functional  relationship  between 
wavelet  coefficient  threshold  values  and  the  desired  compression  ratio  for  several  compression  ratios. 
Other  criteria  for  image  compression  include  optimization  based  on  a  desired  MSE,  or  distortion 
percentage,  will  also  be  investigated.  The  hypothesized  equation  for  image  compression  will  be  re¬ 
evaluated,  with  the  goal  of  minimizing  the  variability  of  any  constants  in  the  equation.  The  sensitivity  of 
the  threshold  selection  equation  to  the  formulation  of  the  wavelet  function  is  also  a  valuable  topic  for 
research. 
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Figure  7.  Image  reconstructions  from  50:1  compressed  representations  generated  by  Argast  (a)'", 
DeVore  (b),  stochastic  optimization  (c),  and  Eq.  21  (d)  (at  56:1)  approaches.  Original  image  (e). 


